Bose-Einstein correlations in WW events at LEP 
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The current status of the LEP results on Bose-Einstein correlations is discussed. Emphasis is 
given to the measurement of Bose-Einstein correlations between decay products from different 
W's, in an energy range between 172 and 209 GeV, dependent on the experiment. For the 
first time all four LEP experiments conclude that no evidence for correlations between pions 
from different W's is seen at the current level of precision. 



1 Introduction 

Correlations between pairs/multiplets of identical bosons (in the simplified experimental prac- 
tice, all like-sign particles are considered instead) are a well known phenomenon, yet the under- 
standing of the effect is far from complete. Let us start with the original observation of the effect 
by Goldhaber and collaborators In order to give an interpretation to their observations the 
authors started from the assumption that contributions from different bosons to the measured 
intensity add incoherently. A strict analogy with the Hanbury-Brown-Twiss effect @ in astron- 
omy was found. In both cases the measured intensity interference reflects the geometry of the 
emitter. However, using this scenario too strictly one soon runs into interpretational problems. 
The shape of the correlation functions measured since Goldhaber's observation until the LEP 
measurements does not reflect the size of the freeze-out volume at the moment of hadronisation 
which extends up to several fermi. Most observations, except in heavy ion collisions, indicate 
source sizes of the order of one fermi. 

An alternative model was proposed by Andersson and RingnerH, in which the correlations 
appear as a coherent effect related to the symmetrisation of the quantum-mechanical amplitude 
corresponding to the full process of particle production in the fragmentation of the Lund string. 
The strong point here is that the introduction of the Bose-Einstein effect becomes less arbitrary 
than before. It essentially depends on two fundamental parameters of the Lund model, the 



string tension k and the hadronisation cutoff b parameter. In this case one can obtain source 
sizes compatible with experimental observations. This model has however one fundamental 
restriction: only bosons from the same string can be subjected to the Bose-Einstein effect, 
provided that there is no Colour Reconnection at parton level. 

For a simple hadronic system like qq from a Z° decay, it may be impossible to decide between 
the two possibilities, since the incoherent approach leaves a freedom of the choice of the input 
particle density, which can be adjusted to reproduce the observed data. 

The study of correlations between two close hadronic systems, such as hadronically decaying 
pairs of WW/ZZ bosons, can eventually help to distinguish between the two possibilities. In 
the incoherent scenario, the difference between correlations within a single hadronic system, and 
correlations between the two systems, should depend only on the overlap of the two systems 
(sources). In the coherent scenario, the correlations between the two systems may not exist at 
all, even for overlapping sources (as long as there is no interaction -colour flow- between these). 

The measurement of inter-W correlations is also important for the estimate of the systematic 
bias in the measurement of the W mass via the direct reconstruction of measured decay products. 
A better understanding of the physical origin of the observed correlations is however necessary 
to ensure a reliable prediction for the uncertainty on the W mass measurement. 



2 Analysis methods 

It is common practice to investigate BEC between particles coming from different W's by means 
of a two-particle correlation function in terms of the Lorenz-invariant four-momentum transfer 

Q = vMpi -P2) 2 : 

where po(Q) represents the two-particle density without the Bose-Einstein effect. This density is 
non-existent in nature and is known as the so-called normalization or reference sample problem. 
We will see that many experiments address this problem in different ways, each with their own 
degree of model and detector dependence. A widely used implementation of the Bose-Einstein 
correlation effect in Monte Carlo generators is the LUBOEI § code, included in JETSET i. 
Experiments use different versions of this code and tune the Monte Carlo samples to their Z° 
data. 



2.1 OPAL analysis 

The OPAL collaboration has published an analysis! for a total collected statistics of 250 pb" 1 . 
In this analysis the two-particle correlation function is constructed using unlike-sign pairs as a 
reference sample and making a double ratio with the correlation function obtained for a Monte- 
Carlo sample without Bose-Einstein correlations at all: 

N data {Q) / N MC {Q) 

W) Nf a l a (Q) / N^(QY { > 

This is done for three samples: fully hadronic WW decays, semi-leptonic WW decays and qq 
events selected as fully-hadronic WW events. One can assume that each of these 3 correlation 
functions can be written as the sum of 3 independent and more interesting correlation functions. 
For example one can write the correlation function for the fully hadronic sample as 

c had (Q) = PLAQ)c s (Q) + Phad(Q) c bJ (Q) + (1 - PfUQ) ~ P^ d (Q))c d (Q), (3) 



where C S (Q), C Z J(Q), C d (Q) represent the correlation functions for particle pairs originating 
from the same W, the Z° background and for pairs originating from different W's, each with 
their own probabilities P(Q), obtained from MC samples without Bose-Einstein Correlations. 
In a next step OPAL makes a simultaneous fit to the three measured correlation function using 
the expression: 

C s > d > z (Q) = N(l + U(Q)X sAZ e- R2Q2 ), (4) 

where /^(Q) is the probability that a given particle pair is indeed a pair of pions, obtained from 
Monte Carlo. Taking into account the distance between the W decay vertices one can impose a 
constraint on the radii: 

(R d ) 2 = {R s f + (correction) 2 . (5) 

This gives a fit result of 

A s = 0.69 ± 0.12(stat) ± 0.06(syst), 
X d = 0.05 ± 0.67(stat) ± 0.35(syst), 

leading to the conclusion that with this method and at the current level of precision it is impos- 
sible to establish whether BEC between different W's exists or not. 



2.2 ALEPH analysis 

The ALEPH collaboration has published results i for the energy range between 172 and 189 
GeV. An update including energies up to 202 GeV was submitted to ICHEP2000 i. Similar 
to OPAL, ALEPH also uses unlike-sign pairs as a reference sample and corrects for resonance 
decays and detector effects by making a double ratio with a MC sample without BEC at all. 
Since the qq background might fake a possible inter-W BEC signal it was decided to add the 
background fraction to the MC reference without BEC included. In this way the two-particle 
correlation function becomes 

r>*( n \ _ N ±± a (Q) , N^ WW+qS) (Q) 

U N data(Qy N MC(WW+qq) (Q) 1 ' 

The distribution of R*(Q) is compared between data and two Bose-Einstein models based on the 
LUBOEI BE3 algorithm, tuned on Z° data, as can be seen in Fig. |]. Fits to this distribution 
for data and models are made using expression 

R*(Q) = K(l + eQ){l + \e- CT2Q2 ). (7) 

The results of the fits are compared by integrating over the correlation signal 

oo 

Xe^ 2 dQ = 4-- (8) 
2 a 

For this measurement ALEPH finds that the value of 1 for the data is compatible with the 
value of I for the BE3 model in which only intra-W BEC are present. The BE3 model with 
intra+inter BEC is disfavored at the level of 2.2 a. 

In a second method mixed semi-leptonic events are used as reference sample. Again a double 
ratio with a MC sample without BEC including the qq background is used, and the two-particle 
correlation function becomes: 

R m (n \ = ly ±± W) / iv ±± WJ_ (q) 
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Figure 1: The R*(Q) distribution for data compared 
with BE3 model predictions. Only statistical errors 
are shown. The solid curve shows the fit result to 
the data. 
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Figure 2: The R m (Q) distribution for data compared 
with model predictions. The solid curve is the fit 
result to the data. 



This distribution (see Fig. ^) is again fitted with a gaussian parametrisation and integrals are 
compared. In this case the inter+intra BEC scenario is disfavored at the level of 3.1 a (stat 
only). 



2.3 L3 analysis 

The L3 analysis relies on a rigorous mathematical treatment and generalizations thereof 111 
and is published □ for the collected data at 189 GeV. A new update has been given for this 
conference, including the 192-202 GeV data. In their formalism one can write the two-particle 
densities for independently decaying W's as 

p ww (l,2) = 2p w (l,2)+2p%%, (10) 



where the second term on the right-hand side of Eq. 10 is obtained by mixing 2 semi-leptonic 
events. In the absence of inter-W BEC the ratio of the left-hand side and right-hand side 
of Eq. 10, which is called D, should be compatible with one. After subtracting 18.6% qq 



background from the fully hadronic term p (1,2), using the LUBOEI BE0 model, L3 makes 
a double ratio by dividing the D distribution for the data by the same distribution obtained 
with a MC sample without any BEC included. This variable is called D' and is fitted with a 
gaussian expression. Both distributions are shown in Fig. ||. The fitted value for the correlation 
strength A is compatible with zero: 

A = 0.013 ± 0M8(stat) ± 0.015(syst). 

Comparison with the inter+intra BEC BE32 model tuned at the Z° data gives a deviation from 
the data of 4.7 a. 



24 DELPHI analysis 

The DELPHI analysis has been updated for a total collected statistics of 531 pb _1 including 
energies from 189-209 GeV. DELPHI used the same formalism as L3 and also studies the 
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Figure 3: The D and D' distribution for data com- 
pared with BE32 model predictions. Only statisti- 
cal errors are shown. Bin-to-bin correlations are not 
considered 
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Figure 4: The D and Ap distribution for the DEL- 
PHI data compared with BE32 model predictions. 
Only statistical errors are shown. Bin-to-bin corre- 
lations are not considered. 



difference between the right-hand side and left-hand side of Eq. 1C, which is called Ap(Q). 
This distribution, together with the D distribution, is shown in Fig. |j. The qq background 
contamination in the data amounts to 15% and is subtracted using the BE32 model tuned on 
Z° data. Since this background subtraction is a delicate point, DELPHI tried to investigate how 
well the BE32 model describes the qq events which are selected as WW events. This was done 
using 4 jet Z° events and high energy qq events with an anti-WW tag. The largest disagreement 
between the model and data did not exceed the 10% level. This study is still ongoing. In order 
to stay as model independent as possible DELPHI does not construct a D' distribution and 
makes a fit directly to the D variable with the following expression: 

D(Q) = N{l + eQ){l + ke- (7Q ) (11) 

After fixing a to 1.01 fm, as was fitted for the inter+intra BEC model prediction, DELPHI finds 
a value of A dat compatible with zero. 

A dat = _ 038 ± 0.057(stat) ± OM(syst) 

The systematic error is still under study and contains for the moment only the contributions 
from the background subtraction (0.05) and from the mixing method (0.03). However, it is 
assumed that these two contributions are the dominant ones. When comparing the fitted value 
of A of the data with the inter+intra BEC BE32 model prediction, DELPHI disfavors the model 
at the level of 3.2 a. 



3 Summary 



It is important to note that for the first time the 4 LEP experiments obtain consistent conclu- 
sions. The LUBOEI models tuned on the Z° data from each experiment, and which include 



BEC between different W's, are excluded by all experiments with varying significance. The LEP 
experiments are on the way to converge on measurement techniques as proposed in EH , which 
is very promising. It is my question to the W-mass measurement community whether they will 
still use these models to estimate their systematic errors. What is clear for me is that WW 
events will not tell us much more about the ongoing discussion on incoherence and coherence, 
and the easy but rather restrictive variable Q might not be the ideal one to be used. Certainly 
a study of multi-string events from LEP1 would be very interesting to address this problem Ej. 
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